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ABSTRACT 

We describe Chandra/ ACIS-I observations of the massive ~ 13-14 Myr-old 
cluster, h Persei, part of the famous Double Cluster (h and x Persei) in Perseus. 
Combining the list of Chandra-detected sources with new optical/IR photometry 
and optical spectroscopy reveals ~ 165 X-ray bright stars with V < 23. Roughly 
142 have optical magnitudes and colors consistent with cluster membership. The 
observed distribution of L^; peaks at L^; ~ xo^°'^ ergs s~^ and likely traces the 
bright edge of a far larger population of ~ 0.4-2 Mq X-ray active stars. From a 
short list of X-ray active stars with IRAC 8 /xm excess from warm, terrestrial- zone 
dust, we derive a maximum X-ray flux incident on forming terrestrial planets. 
Although there is no correlation between X-ray activity and IRAC excess, the 
fractional X-ray luminosity correlates with optical colors and spectral type. By 
comparing the distribution of L^/L^, vs. spectral type and V-I in h Per with 
results for other 1-100 Myr-old clusters, we show that stars slightly more massive 
than the Sun (> 1.5 M0) fall out of X-ray saturation by ~ 10-15 Myr. Changes 
in stellar structure for > 1.5 Mq stars likely play an important role in this decline 
of X-ray emission. 
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Introduction 



Intermediate and low-mass main sequence stars like the Sun show evidence for chro- 
mospheric/coronal activity. Diagnostics of this activity - Ca II H and K emission, X-ray 
emission, etc. - are likely produced from a self-sustaining magnetic dynamo, which results 



from a combination of conv ective energy transport and differential rotation (jParkerl Il955 



Babcocklll96ll : lBarnesll2003l ). Once stars begin to contract onto the main sequence, the X- 



ray luminosity (L^ ) correlates with stellar rotation (v), with L^. oc v^ (jPallavicini et al.lll981 



Giidel et al.l 119971 ). This evolution provides a quantifiable link be tween a major diagnostic 
of activity and the source of activity (see also iPizzolato et al.ll2003l ). Therefore, studying the 
evolution of X-ray emission provides an insight into stellar structure and evolution and the 
influence of activity on the circumstellar environment. 

Young analogs of main sequence field stars have much higher levels of X-ray emission 
than main sequence stars. X-ray surveys of young, < 100 Myr-old open clusters reveal that 
intermediate mass (~ 1-3 Mq) and low-mass stars are strong X-ray emitters. Typical X-ray 
luminosities for young stars are 'y 100-1,000 times larger than the characteristic solar X-ray 
luminosity (IFeigelson et al.ll2005uPreibisch and Feigelsorul2005l : iPreibisch et al.ll2005l ). 



X-ray observations of clusters (e.g. with Chandra) provide detailed constraints on the 
time history of chromospheric activity. As stars age and evolve onto th e main sequence 



their magnetic activity diminishes. Thus, their X-ray luminosity declines (ISkumanichlll972 



Micela et al.l 



1985 



Glide 



Preibisch and Feigelson 



2005f ) 



et al.lll997l : IFeigelson and Montmerlelll999l : lEvans and Sewardll2000 



This decay is likely due to the rotational spin-down of the 
star and reduction of differential rota tion at the radiative/convective zone interface (e.g. 
Noyes et al.lll984l : iBaliunas et al.lll995l ). The youngest (~ 1 Myr old) stars less massive than 
~ 1-2 M0 typically have L^. ~ 10^^-10'^^ ergs s~^ and fractional luminosities (L^ /L^, where 
is the bolometric luminosity) of ~ 1 0~^-10~'^ (IPreibisch and Feigelson! l2005l ). Observa- 



tions of older clusters like the Pleiades ( Micela et al. 19991 ) suggest that the typical X-ray 
luminosities drop to 10^^ ergs s~^ and fractional luminosities drop to ~ 10~^-10~^'^ by 100 
Myr. 

Observations of young stars also show that a stellar mass dependence of X-ray emission 
emerges sometime after ^ 2 Myr. In the COUP observations of the Orion Nebula Cluster, 
young stars have a wide range of fractional luminosities Av hich are not correlated with stellar 
mass ( Preibisch and Feigelson 2005 : Preibisch et al.|[2005l ). X-ray activity from the magnetic 
dynamo is uncorrelated with rotation rate and is described as being " satura ted" , L^./Lj, ~ 
10"^, or "supersaturated", L^,./!^ > 10"^ (IRandich et al.lll996l : ICiidell 12004 ). In contrast, 
observations of older clusters l ike the Pleiades suggest clear relations between L^/L^ and 
stellar mass (IMicela et al.lll999l ). where lower-mass stars have a higher fractional luminosity. 
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Thus, sometime between 1 Myr and 100 Myr, stellar rotation rates decrease su fficiently to 
bring the X-ray activity out of saturation (jPatten and Simon 1996 : Giidel 2004). 



As stars evolve to the main sequence, their X-ray emission may drive important processes 
in the circumstellar disk material from which planets form. X-ray active stars with ages of ~ 
1-10 M yr are surrounded by primord ial circumstellar disks comprised of gas and small dust 
grains (iKenyon and Hartmandll995l ). X-ray irradiation may provide an ionization source for 
primordial disks that po wers disk viscosity to sustain disk accretion when MRI turbulence 
is otherwise ineffective (iGlassgold. Najita. and Igeal 119971 ). Flash heating from X-rays in 
powerful solar fl ares may account for the formation of chondrules in the early solar nebula 
( IShu et al.l 119961 ). X-ray heating of solids during flaring events, perhaps even during the 
early debris disk stage, may also drive chemical rea ctions from spallation which produces 
the short-lived isotopes observed in some meteorites (IFeigelson et al.ll2002l ). 



X-ray irradiation also has important consequences for the evolution of newly-formed 
planets. Large X-ray fluxes can ablate the hydrogen-rich atmospheres of short-period giant 
plan ets (Lammer et al. 120031) and remove the atmosphere entirely for sufficiently high L^. (see 



also lBaraffe et al.ll2004l ). For an X-ray luminosity function scaled to the Pleiades, nearly all 
(~ 85%) Nept une-mass gaseous planets at 0.02 AU from a solar-type star may be eroded to 
Super-Earths (jPenz et al.ll2008l ). G-type stars with high X-ray luminosities can evaporate 
the atmospheres of a wide range of planet masses (~ 0.1-10 M^) within ~ 0.1 AU. In the 
absence of evaporation. X-ray irradiation may also drive a variety o f photochemical reactions 



in the atmospheres of terrestrial planets (e.g. iLammer et al.l 120071 ). 



Understanding the evolution of stellar activity and the X-ray environment of circumstel- 
lar disks and planets requires observations of stars in clusters with ages of ~ 10-30 Myr. Re- 



cent X-ray surveys have concentrated on < 5 Myr-old clusters such as Or ion, NGC 226 
IC 348 and > 50-1 00 Myr old clusters such as a Persei and the Pleiades (iProsser et al. 



1 

Micela et al. 


1999 


2007 


). The large { 



19991 : iPreibisch and Zinneckerl l2002l : iPreibisch and Feigelson 



20051: 



, and 



19961: 



Dahm et al. 



to constrain the evolution of the X-ray luminosity function and the evolution of stellar ac- 
tivity. X-ray data for 10-30 Myr-old clusters also link constraints on the X-ray environment 
during accretion and grain growth in the primordial disk phase with those for > 100 Myr-old 
stars where planet formation should be complete in the terrestrial zone. 

In this paper, we report analysis of Chandra observations of 13-14 Myr old h Persei, 
par t of the famous Double Cluster h and y Persei. Recent Spitzer observations of h Per- 
sei ( ICurrie et al.l l2007al . l2008al : ICurrid l2008bl ) reveal that the cluster harbors a substantial 
population of sta rs with warni dust emission consis tent with debris fr om active terrestrial 
planet formation JCurrie et al.ll2007bl : ICurriell2008bh . At ~ 13-14 Myr Jsiesnick et al.l 12002 



-4- 



Meynet et al.lll993l : ICurridl2008bl ). h Pers ei probes ages during which terrestrial planets are 



poten tially in the final stages of formation (IKenyon and Bromley 



1993') and when nebular gas has recently dissipated (ICurrie et al 



200 



Wetheril 



2007c. see also 



and Stewart 



Thommes. Matsumura. anc 



20081 and references therein), leaving newly-formed gas giants unshielded from high energy 
X-ray photons. 

Our main goals for this paper are twofold. First, we investigate the connection between 
X-ray activity and warm debris emission. For X-ray bright stars with evidence for warm dust, 
we estimate an X-ray fiux in the terrestrial zone that provides input for planetary atmosphere 
evaporation models. Second, we provide constraints on the time-evolution of stellar X-ray 
activity linking studies of younger ~ 1-2 Myr-old clusters like Orion and older clusters like 
the Pleiades. Chandra sensitivity limits likely preclude sampling the full distribution of X- 
ray active stars in h Persei. However, we can analyze the upper envelope of the distribution 
of Lx and L^/L^ as a function of stellar properties. From comparing the bright limit of the 
L^/L^ distribution from h Per with the distribution from younger and older clusters, we look 
for evidence of evolution in stellar activity. A detailed investigation of the X-ray spectra for 
sources with the highest X-ray counts will be included in future work (N. Evans et al., in 
prep.) 



2. Chandra Observations and Ancillary Data 

2.1. Observations, Image Processing, and Source Extraction 

Chandra observations of h Persei were taken with a 41.1 ksec exposure on Dec. 2, 2004 
(Obs. ID 5407; Seque nce Number 200341) with the ACIS detector (chips 0, 1, 2, 3, 6, and 7; 



Weisskopf et al.ll2002l ). The data were obtained in dithered, timed mode, with a frametime 
of 3.2 seconds. On-board event rejection and event telemetry was in VFAINT mode. The 
field was centered on a2nn n=2^19"^00'^, (52ooo=57°07'12", close to the center of h Persei from 



Bragg and KenvonI fl2005f ) (a2ooo=2''18™56.4^ 52ooo=57°08'25") and observed at a roll angle 
of 229°. The data were not registered to an astrometric reference frame (e.g. 2MASS). The 
AClS-1 field covers a 17'xl7' area. 

Figure [1] shows the reduced, smoothed image, binned by 4 pixels and then convolved 
with a 3 pixel gaussian to balance the prominence of the sources in the center of the field 
with those on the edges. On the original data (which is the basis for all further analysis) 
the sources in the center of the field are smaller than sources on the edge where the point 
spread function (psf) is larger. 

We process the images with the standard ANCHORS (An Archive of Chandra Ob- 
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servations of Region s of Star Formation) pipeline (]http : / / cxc . harvard . edu/ ANCHO RS /J 



Spitzbart et al.ll2008l ). Based on the CIAO (Chandra Interactive Analysis of Observations) 
routines, this reduction system is designed to process Chandra images containing many point 
sources and diffuse emission. The CXC pipeline level 2 data products (DS version 7.6) re- 
quired no reprocessing; we only applied an energy filter at 0.3-8.0 keV to eliminate the 
high-energy background. 

To select X-ray sources, we investigated several source detection algorithms available 



in CIAO. We used a recursive blocking scheme with WAVDETECT (iFreeman et al.l 12002 
Mexican hat wavelet detection) with the significance set to yield about 1 false detection per 
field. A 15'xl5' region centered on the aimpoint was extracted at full resolution. WAVDE- 
TECT then identified sources in this region using wavelet scales of 2, 4, 8, 16, and 32. We 
merged the three lists, removing any duplicates. 

This method yields 330 point source detections. Source positions were then adjusted 
using a centroiding routine. We extracted counts using an elliptical extraction region around 
each source based on point-spread function fitting as a function of chip position at 95% 
encircled energy. Background levels were computed using an annular ellipse with the same 
shape as the source extraction. The outer ellipse boundry has major and minor axes six 
times the size of the source; the inner boundry has axes three times the size of the source. 



2.2. Deriving the X-Ray Flux and Temperature 



We derived the X-ray flux and temperature for each Chandra source using a one temper- 
ature APEC model with absorption using the CIAO fitting and modeling package Sherpa, 
V. 3.40. Unless X-ray sources are very strong (100 counts or more), it is useful to fix several 
parameters in the fits to prevent unphysical values. The fixed parameters are the hydrogen 
column density, uh, the initial X-ray temperature, and the abundance parameter. 

We us e the well-cons t rained reddening of h Persei to set n^^. Based on ~ 100 high- 

0.56 ± 0.01, a value w ithin 10% of esti- 



mass stars 



Slesnick et all (120021) derive EfB-V) 



mates from similar work (e.g. iKeller et al.ll200ll : iBragg and KenyonI l2005h and identic al to 
that derived for y P ersei. From optical spectroscopy of ~ 6,000 s tars. ICurrid (l2008bl ) and 



Currie et al.l (j2008d ) derive E(B-V) ~ 0.56 for both h and x Persei. ICurrid ( l2008bl ) identifies 
a full-width half-maximum of ~ 0.1 in the distribution of E(B-V). Very few cluster stars 
have E(B-V) > 0.65 or < 0.45. Because E(B-V) is proportional to hh, adopting a single 



^ http : / / cxc . harvard . edu / sherpaS . 4 / index, html / 
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iiH ~ 3 X 10 cm~ for E(B-V) ~ 0.56 introduces a ~ 10% uncertainty in iih and a ^ 
0.05 dex uncertainty in log(Lx) for each source. Additionally, there is a ~ 10% systematic 
uncertainty as we chose a constant o f proportiona l ity fo r the E(B-V) to N^f rela tion which 
is intermediate between the values of IVuong et al.l (120031 ) and lRyter et al.l (Il996l ). 



We adopted an initial temperatur e kT of 1.5 keV as a compromise between values typical 
of stars i n younger clusters fe.g. M17. iBroos et al.ll2007l ) and stars in older clusters (e.g. the 
Pleiades, iDaniel. Linsky. and Gagnell2002l ). An abundance parameter of 0.3 s olar is routinely 



found in fits of X-ray spectra and assumed constant in our calculations (IFeigelson et al. 



2OO2I ) . The uh and abundance parameters remain fixed throughout the fitting, while kT 
is determined from fitting. The remaining free parameter is the normalization from which 
Sherpa derives both the absorbed and unabsorbed photon fluxes. 

We p erforrned a two-step fitting procedure. First, the unbinned data were fitted using C- 
statistics (jCashlll979l ) and Powell optimization (identifies local fit statistic minimum nearest 
to initial guess. For sources with more than 30 counts, the data were then grouped for 
8 counts/bin and refit using the result from the unbinned fit as an initial guess for the 
normalization and temperature. The binned fitting used x-dvar statistics and Levenberg- 
Marquardt optimization. The details of these optimization methods are described in the 
online Sherpa document at ior0. 

Table [T] lists data for Chandra-detected sources. Column 1 is the source ID, columns 2 
and 3 list J2000 positions, columns 4 and 5 list the raw and net counts in the 95% encircled 
energy ra dius. The har d ness ratios (columns 6-8) describe the energy of X-ray photons. 
Following iGetman et al.l (120051 ). we use three hardness ratios. HRl compares the full range 
of CXO sensitivity: 0.5-2.0 keV vs. 2.0-8.0 keV. HR2 highlights differences in the softer 
region (0.5-1.7 keV vs. 1.7-2.8 keV), while HR3 highlights differences in the harder regions 
(1.7-2.8 keV vs. 2.8-8.0 keV). Columns 10 and 11 list the X-ray temperature and uncertainty 
in temperature. Although a detailed analysis of the X-ray spectra is the subject of future 
work (N. Evans et al., in prep.), the median X-ray temperature for the entire sample is ~ 
1.6 keV, which is reasonable for stars with ages intermediate between ~ 1 Myr and 100 Myr 
as described in §2.2. Columns 12 and 13 list the derived unabsorbed and absorbed X-ray 
fluxes. The median values for the unabsorbed and absorbed fluxes are lO"^"^'^ ergs cm~^ s~^ 
and 10~^^'^ ergs cm~^ s~^. The reduced of the Sherpa fits are listed in column 14. 



^http: / / cxc.harvard.edu/sherpa/mcthods / methods, html 
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2.3. Optical/IR Ancillary Data 



There are several possible types of X-ray active sources on the Chandra field, includ- 
ing chromospherically active pre-main sequence stars, active galactic nuclei (AGN), and 
foreground/background red giant stars. Without useful proper motion data for this distant 



cluster, we rely on optical/infrared colors and deep Chandra number counts (e.g. lBauer et al. 



20041 ) to identify likely cluster members and to make plausibility arguments for sources with- 
out optical IDs. For populous clusters like h and y Perse i, we expect a well-defined locus in an 
optical color- magnitude diagram (e.g. iLyra et al.ll2006l ). The optical luminosities for young 
(e.g 10-50 Myr) intermediate/late spectral type pre-main sequence stars (e.g. G0-M5) are 
up to ~ 5 times greater than their main sequence field star counterparts (e.g. iBaraffe et al. 
19981 ). Thus, field stars have far bluer colors at a given magnitude. Background AGN 



have very blue optical colors compared to pre-main sequence stars and have very faint op- 
tical magnitudes. Foreground red giant stars have far redder colors than cluster stars at a 
given magnitude. The differences in V/V-1 positions for cluster stars and other X-ray active 
sources means that this color-magnitude diagram clearly identifies the positions of X-ra y 
active cluster members in moderately massive clusters such as NGC 2547 (ILyra et al.ll2006l ). 
Because h Persei is > 5-10 times more massive than NGC 2547, the cluster locus should be 
even more easily identifiable. 

To identify the nature of Chandra-detected sources, we merged the Chandra catalog 
with deep VI optical photometry and optical spectra of h Persejf]. Optical photometry of h 
Persei was taken with the Mosaic Imager at the 4-meter Mayall telescope at the Kitt Peak 
National Observator y on October 13-16 and 27-30, 2006, as a part of the MONITOR project 
( lAigrain et al.l 120071 ). Exposures in V and I band were taken using 75 second integrations 
with a 36'x36' field of view cent ered on the cluster. The data were reduce d using the pipeline 
for the INT wide- field survey (jlrwin and Lewis! l200ll : llrwin et al.l 120071 ). correcting for the 
effect s of fringing , cross talk, bias, and atmospheric extinction. Photometry was performed 



as m 



Irwin et al.l (120071 ): instrumental magnitudes were transformed into Johnson- Cousins 
system magnitudes. The catalog contains ~ 42,000 sources detected in at least one band 
and are complete at the 5a level to V ~ 23 and I ~ 20. 

For source matching, we use the optical photometry. For 14 Myr -old cluster stars at 
the distance and re ddening of h Per (E(B-V) ~ 0.56; d ~ 2.3-2.4 kpc ISlesnick et al.l 12002 



Currie et al.ll2008cl ). stars with J-band magnitudes brighter t han the 2MASS completeness 
limit (J ~ 15.7) have V-I < 1.75, V-J < 3.3, and V < 19 JXenvon and HartmannI ^995 



■^The optical pho tometry and spectr oscopy presented in this section will be discussed in more detail in 
an upcoming paper (jCurrie et al.ll2008a) . 
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Baraffe et al. 


1998; 


Siess et al. 


2000) 



data) have V-I ~ 3, V-J ~ 4.5, and J < 18.5 (cf. iKenyon and Hartmannlll995l ). Thus, the 



23 (the 5(7 hmit fo r our optical 



optical data detect fainter, lower-mass cluster stars than the 2MASS data. 

To provide more constraints on the properties of optically-detected Chandra sources, 
we include optical spectra of h Per sources. Optical spectrosco py of ~ 5,000 h Persei sources 



were taken during Fall 2006 and Fall 2007 with Hectospec (IFabricant et al.l |2005| ) on the 



6.5m MMT telescope. We took three 10-15 minute exposures using the 270 g mm~^ grat- 
ing. This configuration yields spectra at 4000-9000 A with 3 A resolution. The data were 
processed using standard Hectospec reduction pipelines. The reduced spectra typically had 
signal-to-noise > 20-40. Finally, we include ~ 1000 archival spectra from the FAST and 
Hydra spect r ograph s on the 1.5 m Till inghast and 3.5 m WIYN telescopes, respectively. 
Currie et al.l (l2007bl ) and ICurrid (l2008bl ) describe these data in detail. The FAST and Hy- 
dra spectra consist mostly of bright, likely early- type stars in h Persei. 

To derive spectral types from these sp ectra, we employed the semi-automatic quantita- 
tive spectral-typing code SPTCLASS from lHernandez et al.l (120041 . see www.astro.lsa.umich.edu/hernandj 
for more information), which is useful for classifying stars regardless of luminosity class and 
surface gravity. SPTCLASS calculates the spectral types of stars using spectral indices 
which compare the line flux of each spectral feature to continuum levels. The relationships 
between spectral indices and spectral types are derived from spectroscopic standards ob- 
served with Hectospec. Thus, SPTCLASS is particularly well suited to spectral type h 
Persei stars. The errors in spectral types are typically ~ 2 subclasses. 

To identify the presence or absence of warm, circumstellar dus t around X-r a y activ e 
stars, we rn a tch AC IS sources to the 2MASS/IRAC catalog from I Currie et al.l (l2007al ). 
Currie et al.l (l2007al ) describe the data reduction, photometry, and source matching between 
2MASS and IRAC in detail. The catalog consists of ~ 31,000 point sources and covers ~ 
0.75 square degrees on the sky, encompassing both h Persei and x Persei. The completeness 
limit for the catalog corresponds to the 2MASS completeness limit (J ~ 15.5); most stars 
brighter than J ~ 15-15.5 have IRAC counterparts in at least one channel (3.6 /xm, 4.5 /xm, 
5.8 /im, or 8 /im). 

Of the 330 X-ray sources, 165 have optical photometric counterparts. The typical po- 
sitional offsets are very small (mean offset is ~ 0.60", median ~ 0.47", and a ~ 0.43") and 
comparable to the astrometric accuracy of Chandra. Of the 165 X-ray/optically-detected 
sources, 101 have spectral types; 123 have also 2MASS/IRAC photometry. We derive the 
X-ray luminosity for optically-detected Chandra sources assuming a distance of 2.4 kpc. 



The catalog of optically-detected Chandra sources along with their net X-ray counts 
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and absorption-corrected X-ray luminositiefl optical/IR photometry, and spectral types is 
listed in Table I 



2.4. Nature of Chandra Sources Lacking Optical Detections 

By comparing the optical survey limits to colors and magnitudes for 14 Myr-old pre-main 
sequence stars and other X-ray active sources we can constrain the properties of Chandra 
sources lacking optical detections. The two main possibilities for these sources are fainter, 
low-mass cluster stars and background AGN. Simple arguments show that both fainter, low- 
mass cluster stars and AGN likely comprise the Chandra-detected population lacking optical 
detections. 

For the Geneva/Bar affe isochrone (see §3.1), stars with V ~ 22 have masses of ~ 0.6 
Mq (Figure [2]). While model uncertainties plague our understanding of the optical col- 
ors/magnitudes of lo wer- mass stars , these stars should be cooler and fainter than 0.6 Mq 
stars of the same age. iBaraffe et al.l (119981 ) predict that h Per stars with ~ 0.25 Mq have 
V ~ 24, two magnitudes fainter than 0.6 Mq stars and undetectable with our optical data. 
If these stars have high fractional x-ray luminosities (> 10~^, Chandra is likely sensitive 
enough to detect them. Because the 5a limit for the optical data is V ~ 23, stars of slightly 
greater mass (< 0.5 Mq) may also lack optical detections but have Chandra detections. 



Assuming a Miller-Scalo Initial Mass Function (IMiller and Scald Il979l ). h Per stars 
with masses between 0.25 Mq and 0.6 Mq are 1.3 times as numerous as cluster stars with 
masses > 0.6 Mq. The number of optically-detected Chandra sources is equal to the number 
lacking optical detections. Therefore, lower-mass cluster stars plausibly comprise a significant 
fraction of the Chandra-detected population lacking optical detections. 



Using the number counts of X-ray sources as a function of limiting flux from lBauer et al. 



( 120041 ). Chandra sources lacking optical detections likely include a significant population of 
optically faint, backgrou nd AGN. The limi ting (unabsorbed) fiux for our survey is log(F^., 
ergs s~^ cm~^) ^ -14.7. iBauer et al.l (120041 ) derive the number of X-ray sources per sq deg 
with fiuxes greater than log (F^. ergs s~^ cm~^) = -14.7 to be roughly 500 soft X-ray sources 
and 1000 hard X-ray sources. Because the ACIS-I coverage is ~ 0.08 square degrees, our h 
Per field should have ~ 40 (80) s oft (hard) AGN. M ost AGN with log F^ ~ -14.7 have faint 
optical magnitudes (R > 22-23; iBauer et al.ll2004l ) and are thus fainter than the 5 a limit 



^AU X-ray luminosities discussed in the following sections are absorption-corrected luminosities 
^The full versions of the tables are available in the electronic edition of this paper 
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of our optical survey. 

Other properties of sources with and without optical Several X-ray sources on the field 
have 2MASS detections but no detectable optical counterparts. These are typically very 
faint and do not define a locus in near-IR color-magnitude diagrams consistent with cluster 
membership. These sources are then more likely to be AGN. 

In summary, the Chandra sources lacking optical counterparts are likely to be a mix 
of low-mass cluster stars and background AGN. Deeper optical data will likely identify the 
origin of some currently unmatched sources. Deeper X-ray data will yield enough counts 
to fit a spectrum to many sources; analyzing their X-ray spectra can distinguish between 
pre-main sequence stars and AGN. 



3. Analysis 



3.1. Optical Properties of Chandra-detected Sources 

Most of the 165 optically-detected X-ray sources define a clear sequence in the V/V-I 
color-magnitude diagram consistent with cluster membership. Figure [2] shows their positions 
(blue and green dots) with respect to all sources (g rey dots) detected within 10' of the h 



Persei center (^2000 ~ 2'^19'"0*, ^2000 ~ 
sources define a sequence from V,V-I ^ 



57°8'35"; see 
' 10,0 to V,V-I 



Bragg and Kenyonll2005l ). The cluster 



24, 3.6, which is clearly separable 
from the background field star population (V > 15 + 3(V-I) for V-I > 1). Foreground stars 
are less numerous than either the cluster stars or the background. The brightest X-ray active 
star (V ~ 8) is located beyond the plot limits. 

To compare the magnitudes and colors of these stars with theoretical predictions, we 
overplot a best-fit isochrone (14 Myr), covering a range in masses from 15 Mq to 0.6 Mq. 
We construct our isochrone for intermediate /low- r nass stars (5-1.4 M m) by interpolating 



between the 12.6 Myr and 15.8 Myr i sochrones fromlBaraffe et al.l (119981). Fo r more massive 



sta rs, we use the Geneva isochrones (iSchaller et al.lll992l : iMeynet et al.lll993l ). On Figure [2l 
the iBaraffe et al.l and Geneva predictions meet at V ~ 17, V-I ~ 1.2. We adjust the isochrone 



to the d istance and reddenin g of h Persei, assuming a distance modulus of 11.85 and E(B-V) 
~ 0.56 (ISlesnick et al.ll2002l ). For the rest of the paper, we label the isochrone constructed 
from tracks for high- mass stars from the Geneva group and intermediate/low-mass stars from 
Baraffe as the ' Geneva/Bar affe' isochrone. 

The agreement between the predicted 14 Myr isochrone and the color-magnitude di- 
agram for all h Per stars is exceptional. The isochrone accurately reproduces the cluster 
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locus from V, V-I = 10, to V, V-I = 22.5, 2.8. At 14 Myr, this range corresponds to stars 
with masses > 0.6 Mq. Isochrones for 12.6 Myr and 15.8 Myr (not shown) also show fair 
agreement, though slightly younger (11.2 Myr) and older (17.8 Myr) isochrones provide a 
much poorer filo. 

The vast majority (142/165) of Chandra-detected sources, especially those fainter than 
V=16, track this sequence extremely well. Many X-ray sources follow the isochrone from V 
~ 10 to V ~ 21. A smaller number of much fainter X-ray sources (V ~ 22-23, V-I ~ 2.75-4) 
also track the locus of cluster stars. In addition to stars in the cluster, 12 foreground stars 
are X-ray active. 

The X-ray active stars also follow clear sequences in an HR diagram (Figure [3l top 
panel) and in a diagram of V-I vs. spectral type (Figure [31 bottom panel). Aside from 
two early-type stars, there is a clear relation between color and spectral type in the lower 
panel of Figure [31 The X-ray active cluster sequence extends from V = 11, V-I= 0.4 to V = 
23.5, V-I= 3.5. Over this range, h Persei sources have spectral types betweeen B2 and MO. 
Most detected X-ray active stars are F5-G5 stars with V=16-20 and V-I = 1-1.5. Early A 
stars with an intrinsic V-I ~ have an observed V-I ~ 0.7 consistent with the mean cluster 



reddening of E(B-V) ^ 0.56 fiBessel and Brettlll988l : ISlesnick et al.ll2002[ ) 



Using the same upper and lower bounds in V/V-I to identify probable cluster members 
lacking X-ray detections, Chandra detects less than 10% of the h Persei cluster members. 
Within 8.5' of the center of the Chandra/ACIS coverage, 1,621 stars with V=ll-22 lie within 
0.75 mags of the isochrone. Of these stars, 130 have Chandra detections. These stars thus 
have optical colors and magnitudes consistent with cluster membership. The total fraction 
of X-ray detected cluster members with the same V, V-I limits and located within 8.5' of h 
Per is ~ 8% (130/1621)0. In comparison, there are 1,764 stars within 8.5' of h Per that lie 
off the isochrone; 28 of these have Chandra detections. Therefore, ~ 1.6% of nonmembers 
show evidence for X-ray activity (28/1764). 

Converting from optical color and magnitude to stellar mass using the Geneva/Baraffe 
isochrone further constrains the fraction of cluster stars detected as a function of stellar 
mass. Very few high mass stars (M=2 M0-5 M©, V/V-I = 0.52 — 0.75) are detected (2.6%, 
6/230). A slightly larger fraction of subsolar-mass stars (0.6 Mq-1 Mq, V-I =1.87 — 2.66) 
is detected (5.9%, 24/407). The detection rate for intermediate (1-2 M©, V/V-I = 0.78 - 



more detailed analysis of the success/failure of is ochrones of different ages as well as isochrones from 
different groups will be included in an upcoming paper (jCurrie et al.ll2008cl ) 

^Given that the locus of giant stars crosses the cluster locus at V-I ~ 1.7, we expect a small amount of 
contamination and thus some uncertainty in the true fraction of X-ray detected cluster members. 
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1.87) mass stars is highest, ^ 11% (97/911). 



3.2. X-ray Luminosity and Fractional X-ray Luminosity 

The vast majority of Chandra-detected sources with optical counterparts have optical 
magnitudes, colors, and spectra consistent with membership in h Persei. We now analyze 
the X-ray properties of likely h Persei members. X-ray active members are defined as those 
within 0.75 magnitudes of the 14 Myr isochron^. We assume an extinctio n of E(B-V) 



0.56, an age of 14 M yr, and a distance modulus of dM=11.85 for the cluster ( ISlesnick et al. 



20021 : ICurriell2008bf ). 



From these membership criteria, we investigate the X-ray luminosity function and the 
distribution of fractional X-ray luminosities. The X-ray luminosity function probes the abso- 
lute X-ray luminosity produced by the stellar chromosphere; the fractional X-ray luminosity 
function explores how X-ray luminous the star is compared to its photosphere. We compare 
these luminosity functions with the stars' dereddened colors and spectral types to investigate 
how they correlate with stellar properties. 

The fractional luminosity, L^./L^, is estimated from the observed L^, and the bolometric 
luminosity, L^,, derived from the well-constrained age, reddening, and distance for h Persei. 
Assuming that all h Persei stars are ~ 14 Myr old and have E(B-V) ~ 0.56, we evaluated 
the expected spectral type, mass, and bolometric luminosity from V,V-I. The derived for 
sources with V-I colors in between points on the Geneva/Bar affe isochrone grid is determined 
by interpolating between the points in theoretical log(L^) space. 



3.2.1. X-Ray and Fractional X-Ray Luminosity Functions 

Figure m shows the observed X-ray luminosity function for h Persei. The number counts 
peak at log(L2;, ergs s~^) ~ 30.3 (2 x 10^° ergs s~^). To compare the p eak to predictions 
for the limiting X-ray luminosity, we follow iFeigelson et al.l ( 20021 . l2005l ) who compute the 



limiting L^. for a given exposure time, source distance, and n^/. For the full 0.5-8 keV band, 
log(Lj. /jm, ergs s"^) ~ 30.3, in agreement with the observed peak. 

We also show the distribution for the subset of stars with V-l=l. 585-2. 225 (dotted line). 



^Technically, the upper limit should be larger than the lower limit because of binarity. For simplicity we 
make the upper and lower limits the same. Using a more restrictive lower limit (e.g. 0.3 magnitudes) does 
not affect our results. 



- 13 - 



At 14 Myr, these stars should have ~ 0.9-1.2 Mq and have G spectral types on the main 
sequence. The distribution for 0.9-1.2 Mq stars is similar to the function for the entire 
population, though it lacks any stars with log(L^.) > 31. 

Because Chandra detects < 10% of all cluster stars, members without detections have 
lower X-ray luminosities. X-ray luminosity functi ons for both young and older clusters hav e 



many stars with log(La; ergs s^^) << 30.3 (e.g. iPreibisch et al.l l2005l : I Jeffries et al.ll2006l ). 
Thus, the distribution of L^. shown in Figure H] identifies the completeness limit in L^., and 
the true X-ray luminosity function peaks at log(La. ergs s~^) << 30.3. 

The distribution of L^/L^ is peaked at L^/L^ ~ 10~'^'^ (Figure |4)d). Because the vast 
majority of cluster stars are undetected, the true fractional X-ray luminosity function is 
likely peaked towards lower values. Three stars have a much higher fractional luminosity 
10~^), and several have a lower fractional luminosity of ~ 10~^'^-10~^ '^. 

Five sources- all early-type, high-mass stars- have an extremely low fractional lumi- 
nosity of ^ lO-'^-lO-^ High -mass stars produce X-ray emission from wind shocks, not a 



magn etic dynamo as in solar /subsolar- mass stars, which yields far lower L^/L^ (jStelzer et al. 



20051 ) . While an unseen low- mass companion could be responsible for the X-ray emission, 
this is an unlikely source for X-ray emission around high-mass stars in h Per. B5 stars 
with L^/U ~ 10^^ typically have V ^ 12; G stars with L^/L^ ^ lO'^-^ have V ^ 18. The 
bolometric corrections for B5 stars are BC ~ -1.6; the corrections for GO stars are BC ~ 
-0.2. Therefore, B5 stars are ~ 7.4 magnitudes brighter than G stars. Thus, the fractional 
X-ray luminosity of a GO star in a G0/B5 binary needed to yield an observed L^/L^ for the 
system of ~ 10~^ is ~ 10^^. Chandra detects only the most X-ray luminous G stars in h 
Per, which have Lx/L^, < 10~^. Because > 90% of cluster G stars are less X-ray luminous, 
the likelihood that an unseen companion generates the X-ray flux observed from B stars is 
low. 

Deeper X-ray observations would yield a better measure of the X-ray luminosity function 
and fractional X-ray lu n iinosi ty function in h Per. In the younger Orion Nebula Cluster, 



Preibisch and FeigelsonI ((20051) detect > 90% of the 0.9-1.2 Mq stars with typical L^/U ~ 
10~^ to 10~^ and a median L^-ZL^ ^ 10~^ to 10~'^'^. At the distance of h Persei, many of 
these stars would lie below the Chandra detection limit for G-type stars. Thus, the true 
X-ray and fractional X-ray luminosity functions are likely peaked at lower and Lx/L^ 
values. 
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3.2.2. X-Ray Properties vs. Other Stellar Properties 

The fractional X-ray luminosity correlates well with dereddened V-I color (Figure [5^). 
Most sources occupy a region bounded by V-I = -0.2-1.3 and log(L^/L^,) ~ -4 to -3. Bright 
cluster stars with blue V-I colors have weak X-ray emission compared to their photospheres; 
fainter stars with red V-I colors have relatively strong X-ray emission. Figures and b 
show that the five extremely low fractional luminosity sources identified in Figure Hb are 
mostly B0-B5 stars. The lack of low-luminosity late-type stars is likely due to the Chandra 
sensitivity limits, which are shown as a dashed line. However, high-luminosity, early-type 
stars (log (Lj./L^) ~ -3) with V-I ~ -0.2 - 0.3 are absent in h Persei and should have been 
easily detected. 

Because brighter clusters stars at a given age have earlier spectral types than fainter 
cluster stars, the fractional X-ray luminosity of h Persei stars also correlates with spectral 
type (Figure Cluster stars with V-I ~ -0.2 - 0.3 have spectral types between ~ BO and 
F5. Though h Persei contains thousands of stars with these spectral types, none of them 
have L,./L, > lO'^ ^ 

In summary. Figure [5] shows that stars with small L^/L^ 10~'^-10~^) are early type 
(B/A), higher-mass stars while those with larger L^j/L^ 10^^-10^^) are intermediate 
type/mass star s (FGK). The B and A stars with X-r ay emission are more massive than ^ 2 



Mq at 14 Myr flBaraffe et al.lll998l : ISiess et al.ll2000l ). The F-K stars with higher fractional 



luminosity have masses of ~ 1-2 Mq. While Chandra likely detects only the most X-ray 
luminous cluster stars, these data provide some constraints on the X-ray luminosity function, 
specifically its upper envelope as a function of spectral type. 



3.3. X-Ray Properties and IR Excess from Warm Debris Disks 



The Do uble Cluster, h and y Persei , contains the largest known population of warm 
debris disks ( ICurrie et al.l l2007a| ]bl. l2008at ICurrid l2008b[). Ongoing terrestrial pla net forma- 
tion is the most likely source for the observed debris (IKenyon and Bromleyl |2004| ) . Here we 
analyze the X-ray properties of Chandra-detected h Persei sources with warm debris disks. 
We focus on a) whether X-ray luminosity correlates with the amount of debris emission and 
b) identifying the X-ray luminosity of X-ray bright stars that are likely forming terrestrial 
planets. 

Figure shows the spectral types of Chandra sources vs. their observed Ks-[8] colors. 
X-ray bright stars with IRAC detections cover a wide range of spectral types (B0-K3). Most 
stars have spectral types earlier than G3 with Ks-[8] ~ 0-0.3. To compare these colors 
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with the predicted colors of stellar photospheres, we add the locus of stellar Ks-[8] colors for 
BO-MO stars from the Kurucz-Lejeune stellar atmosphere models using the SENS-PET tool 
available on the Spitzer Science Center websit^. The locus of photospheric colors is nearly 
independent of spectral type from A stars through G stars (Ks-[8] ~ 0), covering almost the 
entire spectral type range of Chandra sources. 

Figure El also r eveals four Chandra sources with IR excess emission from circumstellar 
dust. ICurrie et al.l (l2007al ) identified 8 /im excess sources as those with Ks-[8] > 0.4 + a [8], 
where a[8] is the typica l uncertainty in the 8 nm magnitude. For FGK stars with typical 
[8] ~ 13.75, a[8] ~ 0.1 (ICurrie et al.ll2007a] ). Three stars have Ks-[8] colors redder than ~ 
0.5 and thus have 8 fim excess consistent with warm, terrestrial zone dust. A fourth star 
without a spectral type also has K5-[8] > 0.5. Thus, at least 4 Chandra detections have 
excess emission at 8 fim. 

The frequency of 8 fim excess emission among Chandra detections is comparable to the 
fraction in the entire cluster. The frequency of Chandra-detected stars with 8 fim excess is 
4/123 (3.3% ± 1.7%. The frequency for Chandra-detected cluster members is 4/106 (3.7% 

.ction ranges from ^ 4% to 8% 



(Currie et al. 


2007a: 


Currie 


2008b: 


Currie et al. 


"1 •'"^ 
2008c) 



bright stars with evidence for warm dust, the similarity in excess frequency between them 
and stars lacking Chandra detections should be considered a tentative result. 

Figure [7] compares the Ks-[8] colors and X-ray luminosities for Chandra-detected cluster 
members. Most stars have log(Lj,., ergs s~^) ~ 30-30.8. Sources with red Ks-[8] colors do 
not identify a unique space in this plot. Three 8 /xm-excess stars have log(La;, ergs s^^) ~ 
30-30.4; a fourth has log(L2;, ergs s~^) ~ 30.8. The Spearman's rank test, a non-parametric 
measure of correlation, reveals that the distribution of with Ks-[8] color is consistent with 
a random distribution (probability ~ 23%, d=0.17). 

Figure [7)d shows the distribution of fractional X-ray luminosity with Ks-[8] color. Most 
IRAC-detected stars have Lx/L^, ~ 10~^ '^^-10~^'^^, which are typical of the X-ray bright 
population as a whole. Sources with 8 nm excess define a slightly more narrow distribution 
between L^./L^ ~ 10~^-10~^'^. According to the Spearman's rank test, the probability that 
Lx/L^, is not correlated with the Ks-[8] color is ~ 0.89% (d=0.44). Thus, there is a 2-3(7 
correlation between X-ray flux and Ks-[8] color. 

If this trend is r eal, the fract i onal X- ray lu minosity probably do es not affect the amount 
of debris emission. ICurrie et al.l (j2007al ) and I Currie et al.l (j2008cl ) find that intermediate- 



'http: / /ssc. spitzer. caltech.edu/tools/senspet/ 
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type (FG) stars more frequently have 8 fim excess than earher stars. Intermediate-type stars 
also have higher fractional X-ray luminosities regardless of their circumstellar environment 
(Figure Eb). Therefore, the correlation between La,/L^ and IR excess in Figure [7|d is likely 
due to the intrinsically larger fractional X-ray luminosities and higher frequencies of warm 
dust for intermediate-type stars compared to early-type stars. 

In spite of these results, we caution that they apply only to the X-ray bright population 
of stars with warm dust, not the population of disk-bearing sources as a whole. Other stars 
with warm dust located within the Chandra coverage must be more X-ray faint. Deeper 
Chandra data are required to more definitively probe the relationship between X-ray activity 
and IR excess. 



3.4. X-ray Active Stars and Planet Formation in h Persei 



If the 8 fim excess emission is due to active terrestrial planet formation, we can estimate 
the X-ray flux incid ent on forming planets given the temperature and location of the dust. 
Currie et al.l (l2008al ) derive dust temperatures of ~ 250-400 K for stars with 8 fim excess. 
Using the bolometric luminosity derived from the Geneva/Baraffe isochrone, we calculate the 
location of the dust assuming that dust grains are in radiative equilibrium with a temperature 
of ~ 300 K. Derived dust locations range from 0.92 AU to 1.36 AU. From the location of 
the dust, we calculate the X-ray flux. 

The incident X-ray flux on objects at the dust locations ranges from 1.53 x 10^ ergs 
cm~^ s~^ to 2.20 X 10^ ergs cm"^ s^^. These levels are equal to or larger than the X-ray flux 
intercepted by the atmosphere of t he hot Jupiter HP 209458b (~ 200 ergs cm"^ s^^ for ^ 
1.1 X 10^^ ergs s~^ at 0.045 AU, see lPenz et al.ll2008l ). Thus, a planet with a mass and radius 
similar to HD 209458b located at ^ 1 AU from an X-ray luminous G star in h Per may 
experience similar rates of atmospheric erosion of up to 
20031 . l2004h . 



10^° g s"^ fiVidal-Madiar et al. 



However, this level of atmospheric esca pe will not l ikely be sustained by any jovian 
planets in the terrestrial zones of h Per stars. iPenz et al.l (120081 ) indicate that the timescale 
for erosion of a signiflcant fraction of a planet's mass is on the order of 100 Myr-1 Gyr. Over 
this timescale, the X-ray lurn inosity of the star will drop by several orders of magnitude 
(jPreibisch and Feigelsonll2005l ). which greatly reduces the rate of atmospheric escape. 
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4. Comparison With Other 1—100 Myr-old Chandra-observed Clusters 



X-ray emiss ion from young stars is most often associat ed with a magnetic dynamo, the 
"a-fi" dynamo (Parker 1955 : Babcock 1961 : BarnesI boosl ) . Differential rotation between 
an inner radiative core and an outer convec tive envelope driv es the dynamo and powers 
chromospheric and coronal activity. For the ISiess et al.l (120001 ) isochrones of a 14 Myr-old 
cluster, we expect stars with outer convective envelopes to have spectral types later than ~ 
GO and stellar masses less than ~ 1.4-1.5 Mq. 

Early-type stars without convective envelopes do not produce a dynamo. Thus, these 
stars do not display stellar activity. However, many early-type stars have large X-ray lu- 



mino sities with typical Lx/L^, ^ 4 — 400 x 10 ^ (e.g., IChlebowski et al.lll989l : IStelzer et al. 



20051 ). Shocks in high velocity stellar w inds are the most likely source of X-ray emission in 



these stars (e.g., iFeldmeier et al.l 119971 ). For h Per, we expect wind emission to dominate 



chromospheric activity for stars with spectral types earlier than FO and masses larger than 

~ 1.7 Mq. 

Our results suggest that h Per contains both types of X-ray sources. Late-type stars 
likely powered by the a-Vt dynamo have Lx/Li, < 10^^ (Figure [5]). Stars with spectral 
types earlier than ~ FO have fractional X-ray luminosities comparable to typical wind-driven 
sources. These fluxes are probably not produced by an unseen late-type companion (§3.2.1). 
Thus, shocks probably produce X-ray emission in these stars. 

In this section, we compare the X-ray properties of h Per stars with X-ray active stars 
in younger and older clusters studied by Chandra and XMM-Newton. We divide our sample 
into "early" (B0-A3) and "intermediate" (FGK) type stars. Comparing h Per data with 
results for other clusters then provides constraints on the evolution of X-ray activity from 
intermediate-type stars and stellar wind-driven emission from early-type stars. 



4.1. X-Ray Emission from Early- Type, High-Mass Stars (BO- A3) 



Stelzer et al.l (120051 ) divide early type X-ray active stars in the ONC into two groups. 



"Strong wind" stars have spectral types between O and B3 (M > 7 Mq); "weak wind" stars 
have spectral types between B5 and A9 (M ~ 2-7 M©). More than half (16) ONC early-type 
stars show clear evidence for X-ray emission. Orion stars with spectral types earlier than B3 
are strong-wind stars with X-ray emission. About 64% of stars with spectral types between 
B5 and A9 are detected by Chandra. 



At 14 Myr, assuming that the division between strong and weak-wind stars is set by stel- 
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lar mass, strong wind sources in h Per have spectral types earlier than B3. Weak wind stars 
have spectral types between B3 and A3. Chandra detects 9 h Per stars ea rlier than A3. None 
of the se stars correspond to Be sta rs or candida t e Be st ars identified by iBragg and Kenyon 
( I2OO5I ) , ISlesnick et al.l (120021 ) , and ICurrie et al.l (l2008al ) and none show Hq, emission char- 
acteristic of Be stars. Compared to the total population of B0-A3 stars within 8.5' of h 
Persei (where V-I < 0.76), ~ 4% are detected by Chandra. About 5.7% (6/106) of stars 
earlier than B3 (V-I < 0.4) are detected, while 2.8% (6/217) of stars between B3 and A3 are 
detected. Other early-type cluster members either lack X-ray activity or have wind-driven 
X-ray luminosities less than ^ 10^^ ergs s~^. 

If the X-ray luminosity function for strong and weak-wind sources does not evolve from 
~ 1 Myr to ~ 14 Myr, Chandra should have detected more weak wind and strong-wind h 
Per stars. Most (6/9) of the strong- wind stars in Orion have X-ray luminosities greater than 
the peak in Lx for h Per (~ 30.3 ergs s^^). Similarly, 44% (3/7) of the weak-wind Orion 
stars have > 30.3 ergs s^^. Although sensitivity limits for h Per observations preclude 
more detailed comparisons between the ONC and h Per luminosities for early-type stars, 
these results may suggest that the X-ray luminosities for early-type, high-mass stars at 14 
Myr are weaker than they are at 1-2 Myr. 



4.2. X-Ray Emission from < 3 M© (FGK) Stars 



Recent Chandra and XMM-Newton observations of young stellar clusters provide some 
constraints on the evolution of X-ray activity as a function of age and stellar mass for 
intermediate and low-mass stars. Independent of stellar mass, the typical fractional X-ray 
luminosity of 1 Myr-old stars in the ONC is L^^/L^ ^ 10~^ (jPreibisch et al.l l2005l ). For 
most young stars, the st ellar rotational velocity incre ases with stellar mass, so Lx/L^ is 
independent of rotation (jPreibisch and Feigelsonl |2005| ) . By analogy with dynamo models 
for X-ray activity, the fractional X-ray luminosity for these stars is " saturated" . 

By 40-50 Myr, when intermediate-mass stars are on t he main sequence, the stellar ro- 
tatio n rates of intermediate- mass have declined significantly ( Prosser et al.lE996: 1 Jeffries et a. 



1996 



20061 ) and correlate better with X-ray fiux. Stars in 50 Myr-old a Persei ( iProsser et al. 
Fig. 5) show a systematic decline in Lx/L^ from 1.2 M© (B-V ~ 0.6, Gl at 50 Myr) to 1.4 
Mq (B-V ~ 0.4, F5 at 50 Myr). Cluster stars with masses < 1.2 M© have Lx/U 



10 



-3 



and thus have saturat ed X-ray emission. The distribution of Lx/L^ for NGC 2547 (38 Myr) 
shows a similar trend (I Jeffries et al.ll2006l ). 



At ages > 100 Myr, subsolar-mass stars leave X-ray saturation. A wide range of stars 
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in the Pleiades (M < 0.85 M©; K3 at 100 Myr) are intrinsically much fainter in X-rays than 
younger stars with the same mass. Finally, older field stars more m assive than > 0.1 M^ti 



also exhibit a decline in X-ray luminosity relative to younger clusters ( iPreibisch et al.l 12005 
Preibisch and Feigelsoru |2005| ). L^/L^ ranges from 10~^-^ to 10~^ for solar-mass stars and 
10"^ to 10"^-^ for subsolar -mass stars. The evolution of fractional X-ray luminos ity vs. time 



show s a systematic decline and is consistent with the ISkumanichI braking track (jSkumanich 
19721 . L^(t) oc t-i/2). 



At least two mechanisms are responsible for the evolution of L^/L* as a star contracts 
onto the main sequence. Rotational spin-down of a star reduces differential rotation at 
the radiative zone/convective zone boundary. Because differential rotation is responsible 
for driving the magnetic dyna mo that produces X-ray emission, X-ray activity diminishes 
( iPreibisch and Feigelsoru l2005l ). 



Intermediate and high-mass stars also undergo significant structural changes as they 
evolve onto the main sequence. At the age of the Orion Nebula Cluster 2 Myr), mai n 
sequence A stars (1.7-2.5 Mq) have spectral types between Kl and K4 (ISiess et al.l 120001 ). 
From 2 Myr to 30 Myr, ~ 2 Mq stars evolve in spectral type from K3 to G8 (5 Myr) to 
A3 (30 Myr). Main sequence F3 stars (~ 1.5 Mq) also undergo substantial changes in their 
internal structure, evolving from at ~ K5 at 2 Myr, to KO at 10 Myr, to F3 at 30 Myr. 
As stars evolve from G/K stars to A/F stars, their outer convective zones shrink and then 
disappear. The disa ppearance of their convective zones eliminates magnetic dynamo- induced 
X-ray activity (see 



Gimiandlll986l : IStelzer et al.ll2005l : IPreibisch and Feigelsonll2005[ ). 



The X-ray detections for stars in h Per place useful constraints on both of these mech- 
anisms. For h Per stars with masses ~ 1.0-1.5 Mq, the maximum in log L^/L^, (~ —2.8 
to —3.0) is similar to the maximum observed in the ONC. Thus, there is no clear evidence 
for evolution in X-ray activity among 1-1.5 Mq stars from 1 Myr to ~ 14 Myr. Deeper 
Chandra data are required to probe the X-ray luminosity function of solar and subsolar-mass 
stars at smaller Lx/L^, and to constrain the evolution of their X-ray activity. 

For more massive stars, there is clear evidence of evolution in X-ray activity. At ^ 14 
Myr, h Per stars with AO to ~ GO spectral types have masses of 3 Mq to 1.5 Mq. In the 
ONC, stars with these masses have a maximum log L^/L^ f» —2.7 to —3.0. This maximum 
is independent of spectral type. In h Per, the maximum L^/L^, clearly declines with spectral 
type. Our Chandra data yield an approximate relation log {L^/ L^)max ~ — 3.5 — (M^/1.5M0) 
for 1.5-3 Mq stars. Thus, the fractional X-ray luminosity declines by a factor of roughly 30 
(3) for 3 (1.5) Mq stars. 



Pre-main sequence stellar evolution is the simplest explanation for the evolution of log 
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[Lx/ Li,)max in 1.5-3 Mq stars. As these stars evolve to the main sequence, their effective 
temperatures increase and they lose their convective envelopes. Stars with masses ^1.8 Mq 
have spectral types earlier than A7 and lack the convective atmospheres need for magnetic 
dynamo-driven X-ray activity. Thus, the X-ray activity of these stars may be similar to the 
weak wind sources at earlier spectral types. Stars with masses ~ 1.5-1.8 Mq have ~ A7-G0 
spectral types and much smaller convective atmospheres than G-type stars. These stars can 
drive weak dynamos and thus have larger fractional X-ray fluxes than earlier type stars. In 
both sets of stars, the decline in X-ray activity appears linked to changes in stellar structure 
as the stars evolve from K-type stars at ~ 1-2 Myr to A-type stars on the main sequence. 

Deeper Chandra data in h and x Pg^' ^-nd other 5-50 Myr-old clusters can test these 
conclusions. If stellar evolution is responsible for the decay of X-ray emission among more 
massive stars approaching the main sequence, this decay should be correlated with spectral 
type and uncorrelated with stellar rotation. Chandra observations of other young clusters 
are needed to see how {L^/ L^)„iax and the X-ray luminosity function of massive stars evolves 
with cluster age. If changes in the magnetic dynamo are responsible for the decay of X-ray 
activity among lower mass stars, deeper observations of (i) 10-20 Myr clusters (to measure 
the X-ray luminosity function for 0.1-1.5 Mq) stars) and (ii) 20-50 Myr clusters (to measure 
evolution in {L^/ Li,)max) are required. 

5. Summciry 

We analyze the first Chandra survey of h Persei, part of the Double Cluster in Perseus. 
By matching Chandra-detected sources with optical/IR photometry and spectra, we identify 
X-ray active cluster members. Using the well constrained age and distance to h Per, we 
compute the bolometric luminosity of h Per stars and compare this luminosity to the derived 
X-ray luminosity. Comparing the X-ray luminosity and fractional X-ray luminosity to optical 
colors, spectra, and near-to-mid IR colors, we investigate the connection between X-ray 
activity and stellar /circumstellar properties. 

These analyses yield the following results: 

• There is a clear correlation of V, V-I with spectral type for the X-ray detections. The 
position of X-ray active stars relative to the observable cluster locus in V, V-I and 
their locus in V vs. spectral type shows that the vast majority of these stars are likely 
cluster members. 

• Within the ACIS coverage, Chandra detects 8% of cluster stars more massive than ^ 
0.6 M„. About 10% of 1-2 Mq stars are detected while ~ 3% of 2-5 M„ stars and 
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~ 6% of 0.6-1 M0 stars are detected. Thus, the sources detected by Chandra hkely 
correspond to the bright end of a much larger population of X-ray active stars. 

The observed distribution of X-ray luminosities for Chandra-detected h Per sources 
peaks at log(L2^) ~ 30.3. Because of the sensitivity limitations of Chandra, the true 
X-ray luminosity function likely peaks at a lower luminosity. 

The observed distribution of fractional X-ray luminosities correlates well with the V-I 
color and with spectral type. Early-type stars (M > 2 M©) have a very small fractional 
luminosity (log(La,/L^) < -4 to -5); L^./!^ is much larger for G and K stars (M ~ 1-1.4 
Mo, (log(L,/Lj ~ -4 to -3). 

There is no evidence for an X-ray saturated population of stars more massive than 
~ 1.5 Mq. In the Orion Nebula Cluster, the most X-ray active 1.5-3 Mq stars have 
log(Lj./Lj,) -3 independent of stellar mass. In h Persei, the most X-ray active 1.5-3 
M0 stars have log(L2./L^) ^ -3.5 - M^/1.5M0. Changes in stellar structure as these 
stars contract onto the main sequence likely play a critical role in this decline of X-ray 
activity. 

From our sample of X-ray bright stars, the presence of warm, terrestrial zone circum- 
stellar dust is not correlated with X-ray luminosity. While stars with higher fractional 
X-ray luminosity are more likely to have warm dust, this trend is likely due to the 
higher frequency of warm dust around intermediate-type stars than around early-type 
stars. 

The warm circumstellar dust emission from X-ray active stars studied i n this paper 



likely res ults fro m debris emission due to terrestrial planet formation (e.g. iCurrie et al. 
2007al lbl. [2008al Hl By computing the X-ray luminosity of these sources and assuming 
a dust temperature of ~ 300 K, we derive an incident X-ray flux on forming terrestrial 
planets in systems with X-ray bright stars. This flux is very small compared to stellar 
irradiation and the radiated flux from terrestrial zone dust. Thus, X-ray emission is 
unlikely to affect the detectability of warm debris disks. 



These data show that h Persei is important in understanding the evolution of X-ray activity 
for a range of stellar masses and provides connections between X-ray activity and terrestrial 
planet formation. Deeper Chandra observations of h Persei can better constrain the X-ray 
luminosity function to lower stellar masses, provide a much larger sample from which to 
compare X-ray activity and planet formation, and more clearly reveal the characteristic X- 
ray environment experienced by forming terrestrial planets. Deeper Chandra observations 
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will also yield enough counts for X-ray bright sources to study the evolution of coronal 
temperatures in intermediate-mass stars. 
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Table 1. Chandra/ACIS Detections 



ID RA DEC Raw Counts Net Counts HRl 


HR2 


HR3 UH 
J022 


kT CT(kT) Flux„, 
cm~^ (keV) (keV) ergs s' 


nahs. 
~^ cm 


— 2 —1 

ergs s cm 


Red.(x^) 

-2 


1 34.9295 TjT.lOOr, 34 


31.7 


-0.1593 


-0.2562 


0.0054 0.3 


3.903 -999 5.49E- 


15 


3.81E-15 


0.572 


2 34.9195 57.0918 6 


3.2 


-999 


-999 


-999 0.3 


0.216 0.044 9.38E- 


15 


1.2BE-15 


0.302 


3 34.9162 67.0697 28 


20.3 


-999 


-0.7477 


-999 0.3 


0.913 0.429 4.25E- 


■15 


1.84E-15 


0.693 


4 34.9039 67.0647 26 


20 


-0.5233 


-0.0469 


-0.6703 0.3 


1.64 0.641 4.01E- 


■15 


2.15E-15 


0.487 


5 34.8843 67.1138 17 


15.2 


-0.7636 


-0.6217 


-0.4442 0.3 


1.699 0.14 2.54E- 


■15 


1.38E-15 


0.632 



Note. — Sources detected by Chandra/ACIS. The hydrogen column, rvjq, is fixed at 3X 10"^^ cm ^ based on the well-constrained reddening to 
h Persei. The temperature (kT), normalization, and x-ray flux {^^x^'X-unahs.) ^'^^ determined by Sherpa. 



Table 2. Catalog of Chandra sources with Optical Counterparts 



ID RA DEC Net Counts (T(Net Counts) log(L^) ergs 3"^ ST o-(ST) ,T H K,, [3.6] [4.5] [,5.8] [8] V tj(V) I <t(I) 



3 


34.9151 


57.0597 


20.3 


5.83 


30.47 


99 


99 


14.026 


13.536 


13.340 


13.144 


13.121 


13.151 


13.425 


16.512 


0.001 


15.050 


0.001 


4 


34.9035 


57.0546 


20.0 


6.74 


30.46 


48.6 


1.4 


15.883 


15.235 


14.932 


14.870 


99 


99 


99 


19.094 


0.004 


17.309 


0.003 


5 


34.8843 


57.1137 


15.2 


4.36 


30.25 


30.5 


2.2 


13.897 


13.521 


13.473 


13.246 


99 


13.425 


13.274 


15.872 


0.001 


14.785 


0.001 


7 


.'i4.870<> 


.-•7.1.-1 4.-1 


11).:', 


2.1)6 


:if).17 


i:i.2 




11.044 


l().<>7.-i 


1().<)7() 


1().!)2() 


If). 1)24 


ILDS.^ 


11). 1)84 


11.X14 


0.001 


11.174 


0.001 


8 


.■i4.87()."i 




17.2 


4.!)4 


:!U.:!!) 


!)!) 


<J<) 


l.-|.(i24 


14.<).-|1 


14..-|<).-| 


14.4.-i2 


14.462 


!)!) 


!)!) 


18.771 


0.00:i 


l(i.982 


0.003 



Note. — Chandra sources with optical counterparts. The spectral types (ST) correspond to numerical spectral types: 10=B0, 20=A0, 30=F0, 
40— G(], and 5(] — KO, etc. The uncertainties in spectral types ((t(ST)) are given in subclasses. Sources without spectral types and/or photometry 
in a given filter have entries denoted with '99'. 
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Fig. 1. — False-color image of the Chandra field binned by 4 pixels and convolved with a 
3- pixel gaussian. The center of h Persei is at 0:2000 ~ 2^18'"56.4*, ^2000 ~ 57°8'25" according 
to iBragg and Kenyoru (120051 ). 
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Fig. 2. — V/V-I color-magnitude diagram with h Persei sources within 10' of the cluster 
center (grey dots). Blue dots correspond to Chandra-detected sources with Ks-[8] < 0.5 
and light-green dots correspond to Chandra-detected sources with Ks-[8] > 0.5. The solid 
line denotes the best-fit (14 Myr) isochrone using the Geneva/Baraffe models with 0.75 
magnitude upper and lower limits. The lower mass limit for the isochrone (V=22, V-I=2.8) 
corresponds to 0.6 Mq. Most Chandra sources with optical counterparts follow the predicted 
isochrone. 
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Fig. 3. — V vs. spectral type (top panel) and V-I vs. spectral type (bottom panel) for 
Chandra-detected sources. The sequence of X-ray active stars clearly traces a locus in 
both diagrams. We overplot error bars for both V-band photometry and spectral type. 
Photometric errors are too small to see for most sources. 



-31 - 



60 



40 - 



20 



30 f 



30.0 



30.5 31.0 
og(Lx) ergs/s 



31.5 



20 - 



10 - 



-5 -4 -3 

og(Lx/Lstar) 



Fig. 4. — Observed distributions of and L^^/L^ for the h Per stars with optical/IR counter- 
parts, (top) Histogram of the for Chandra sources. The number counts peak at log(Lj;) 
= 30.3 ergs s"-*^ for all stars (sohd line) and at log(L3;) ~ 30-30.4 ergs s~-^ for stars with M 
1 Mq (dashed hues) . (bottom) Histogram of the fractional X-ray luminosity. The fractional 
luminosity is peaked at log(L3;/L^) ~ -3.5 for all stars (sohd line) and at -3.2 for stars with M 
~ 1 M0 (dashed hne). X-ray radiation contributes ~ 10~^-10~^ of the total stellar flux for 
most h Per stars. In both plots, Chandra sensitivity limits preclude us from fully sampling 
the entire X-ray luminosity function and fractional X-ray luminosity function. 
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Fig. 5. — Fractional X-ray luminosity as a function of dereddened V-I color (top panel) 
and spectral type (bottom panel). The X-ray luminosity defines a clear sequence from 
bluer stars with a low luminosity to red der stars with a highe r luminosity. The nominal 
Chandra detection limit determined from iFeigelson et al.l (120051 ) is shown as a dashed line. 
Uncertainties for luminosity and fractional luminosity shown in both panels are determined 
from the Sherpa fits; uncertainties for the color and sp ectral type are also sho wn. Stellar 



masse s for 14 Myr-old stars from iBaraffe et al.l (1l998l ) for low-mass stars and ISiess et al. 
(120001 ) for high-mass stars are indicated by vertical dotted lines. 
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Fig. 6. — K5-[8] vs. spectral type for optically-detected Chandra sources consistent with 
cluster membership. The dashed line shows the locus of photospheric Ks-[8] colors from the 
SENS-PET tool (available on the Spitzer Science Center website) using the Kurucz-Lcjcune 
stellar atmosphere models. The dotted line shows the division between photospheric and 
excess sources adopted here. Error bars for excess sources are shown exphcitly; the error 
bar in the lower left corner indicates a typical uncertainty in the K^-fS] color {a ~ 0.1) for 
all sources. Cooler stars have a larger range of -[8] colors. At least three of the X-ray 
active F and G stars have red 2MASS-IRAC colors consistent with warm, terrestrial zone 
dust emission. 
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Fig. 7. — The relationship between X-ray luminosities and IR excess emission. L^, vs. K^- 
[8] (top panel) and log(La;/L*) vs. K5-[8] (bottom panel) for sources detected by Chandra, 
2MASS, and the 8//m filter in IRAC. The dotted fine corresponds to the division between 
excess emission and no excess emission. The uncertainties in (fractional) luminosity and 
photometry are shown for excess sources; the typical Kg- [8] uncertainty is also shown (lower 
left-hand corner). Most sources lie in a distribution between Ks-[8] =0 and 0.4. About 3-4 
stars may have IRAC excess emission indicative of warm, terrestrial zone dust. 



